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ABSTRACT: One of the major challenges in prostate cancer 
therapy remains the development of effective treatments for 
castration-resistant prostate cancer (CRPC), as the underlying 


mechanisms for its progression remain elusive. Previous studies 
showed that androgen receptor (AR) is crucially involved in 
regulation of metabolism in prostate cancer (PCa) cells 
throughout the transition from early stage, androgen-sensitive 
PCa to androgen-independent CRPC. AR achieves such 
metabolic rewiring directively either via its transcriptional 
activity or via interactions with AMP-activated protein kinase 
(AMPK). However, due to the heterogeneous expression and 
activity status of AR in PCa cells, it remains a challenge to 
investigate the links between AR status and metabolic 


alterations. To this end, we compared the proteomes of 


three pairs of androgen-sensitive (AS) and androgen- 
independent (AI) PCa cell lines, namely, PC3-AR*/PC3, 
22Rv1/Dul45, and LNCaP/C42B, using an iTRAQ labeling 
approach. Our results revealed that most of the differentially 
expressed proteins between each pair function in metabolism, 
indicating a metabolic shift between AS and AI cells, as further 
validated by multiple reaction monitoring (MRM)-based 
quantification of nucleotides and relative comparison of fatty 
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acids between these cell lines. Furthermore, increased adenylate kinase isoenzyme 1 (AK1) in AS relative to AI cells may result in 
activation of AMPK, representing a major regulatory factor involved in the observed metabolic shift in PCa cells. 
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E INTRODUCTION 


Prostate cancer (PCa) is the most prevalent cancer among 
males.' PCa patients are initially treated primarily with 
androgen-deprivation therapy (ADT). However, most cases 
of advanced prostate cancer progress to castration-resistant 
prostate cancer (CRPC) after a period of ADT. Due to the 
limited use of chemotherapy in the management of CRPC, 
there remains an urgent need for effective therapy that delays or 
prevents the landmark events that characterize the morbidity 
associated with this cancer.” 
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PCa is very heterogeneous in terms of grade, genetics, ploidy, 
and oncogene/tumor suppressor gene expression, and its 
biological, hormonal, and molecular characteristics are ex- 
tremely complex.** Previous studies on clinical samples and 
xenograft models of PCa indicate that abnormal activation of 
androgen signaling plays a pivotal role in the progression of 
CRPC.° As a result, the sensitivities of different PCa cell lines to 
androgen are heterogeneous. The differing androgen sensitiv- 
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ities of PCa cell lines are well-documented (http: //capcelllines. 
ca/) and represent the multiple stages of PCa progression.’ Cai 
et al.° suggested that, in CRPC, partial restoration of the 
transcriptional activity of the androgen receptor (AR) allowed 
active expression of the genes that mediate androgen 
biosynthesis and maintained intracellular androgen level as 
well as AR activity. Recently, abiraterone acetate, an androgen 
biosynthesis inhibitor, was shown to improve overall survival in 
patients with CRPC after chemotherapy and also showed a 
trend toward improved overall survival in patients who had not 
received previous chemotherapy.” 

Advanced CRPC patients treated with abiraterone acetate 
generally relapse within a year, and AR appears to be active in 
the relapsed tumors.® AR signaling has been shown to cross-talk 
with alternative survival pathways, such as PI3K/AKT and 
MAPK pathways.” This may partially explain the resistance to 
AR-targeted therapies in CRPC. However, interactions between 
AR signaling and these survival pathways are cell-type- and 
context-dependent.'” Hence, conclusions based on one PCa 
cell line could not be directly transferred to all other PCa cell 
lines or applied to clinical situations. 

Here, we set out to obtain a systematic and unbiased 
understanding of the mechanisms underlying CRPC pro- 
gression. To this end, we used the iTRAQ labeling approach 
and performed pairwise comparisons between three androgen- 
sensitive (AS) and three androgen-independent (AI) PCa cell 
lines. Also, multiple reaction monitoring (MRM) analyses were 
applied to targeted metabolites in order to validate the 
metabolic changes assumed. Our results indicate that a 
metabolic shift and distinct forms of cell adhesion existed 
between AS and AI PCa cells. 


mM EXPERIMENTAL SECTION 


Cell Lines and Reagents 


Prostate cancer cell lines PC3, Dul45, C42B, 22Rvl, and 
LNCaP were purchased from the American Type Culture 
Collection (USA). LNCaP cells were maintained in RPMI- 
1640 supplemented with 10% fetal bovine serum (FBS), 100 
U/mL penicillin, and 100 mg/mL streptomycin; the other cells 
were maintained in DMEM supplemented with 10% FBS, 100 
U/mL penicillin, 100 mg/mL streptomycin, and 2 mM 
glutamine. The PC3-AR* cell line, constructed by stable 
expressing of GFP-tagged AR in PC3 cells, was kindly provided 
by Prof. Gao Wei-giang. Cell lines were cultured at 37 °C in a 
humidified atmosphere containing 5.0% CO). 

Forty two methyl ester standards of fatty acids and 24 kinds 
of fatty acids were purchased from Cayman (Michigan, USA) 
and Sigma (Missouri, USA) as listed in Supporting Information 
Table S1. Eight standards for nucleotides and acetyl-CoA were 
purchased from Sigma and are listed in Table S1, and 
commonly used chemical reagents were purchased from 
Sigma. All chromatographic reagents were from Avantor 
Performance Materials (Pennsylvania, USA). Trimethylsilyl 
diazomethane was purchased from Acros Organics (Massachu- 
setts, USA). 


Protein Extraction and Denaturing 


Cells were washed with cold phosphate-buffered saline (PBS) 
twice and harvested by scraping in ice-cold lysis buffer (8 M 
urea, 100 mM Tris-HCl pH 8.5, and protease inhibitor cocktail 
(Roche, Basel, Switzerland)). Cells were sonicated on ice, and 
cell lysates were centrifuged at 16 000g for 15 min at 4 °C to 
remove cell debris. Supernatants were collected and stored at 
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—80 °C until use. Cell lysates were denatured according to the 
FASP procedure with minor modifications.'’ Briefly, the cell 
lysates were transferred to Amicon-0.5 ultrafiltration units, and 
the buffer was exchanged with UA buffer (8 M urea, 0.1 M Tris, 
pH 8.5) three times. Protein samples were reduced with 10 
mM dithiothreitol at 37 °C for 1 h and then alkylated with 50 
mM iodoacetamide at room temperature for 1 h (in the dark). 
The buffer of the denatured protein samples was then 
exchanged with UA buffer three times and exchanged with 50 
mM tetraethylammonium bromide (TEAB) twice. Proteins 
were recovered by reverse centrifugation at 1000g for 2 min and 
resolved in 50 mM TEAB. Protein concentrations were 
determined in triplicate using the bicinchoninic acid method 
(Thermo Scientific Pierce, Massachusetts, USA) according to 
the manufacturer’s instructions. 


Protein Digestion and iTRAQ Labeling 


Equal amounts of protein (320 ug) were digested with trypsin 
at 37 °C overnight. Then, 80 yg of tryptic peptide was labeled 
with 8-plex iTRAQ reagent according to the manufacturer’s 
instructions. The six cell lines were labeled as follows: PC3 with 
113 tag, PC3-AR* with 114 tag, C42B with 115 tag, LNCaP 
with 116 tag, 22Rvl with 117 tag and 119 tag (technical 
replicate), and Dul45 with 118 tag and 121 tag (technical 
replication). Following labeling, equal volumes of labeled 
samples were aliquoted and mixed to test the labeling efficiency 
by LC-MS/MS analysis. 


Peptide Fractionation Using High pH Reversed-Phase 
Liquid Chromatography 


The labeled peptide samples were pooled and desalted using 
Sep-Pak Vac 1 cm? (100 mg) C18 SPE cartridges (Waters, 
Massachusetts, USA). Briefly, the column was washed with 
buffer B (80% ACN, 50 mM NH,HCO,) three times and 1 mL 
each time, then washed with buffer A (50 mM NH,HCO,) 
three times and 1 mL each time. Samples were loaded onto the 
column three times; the column was washed with buffer A three 
times and 1 mL each time, then eluted with 0.5 mL of 40% 
ACN, 50 mM NH,HCO,, and 0.5 mL 60% ACN, 50 mM 
NH,HCO,, sequentially. The combined elute was concentrated 
in a CentriVap benchtop centrifugal vacuum concentrator 
(Labconco, Missouri, USA) at 30 °C. The volume was taken to 
33 uL with solution A used in peptide fractionation. Thirty 
microliters of the sample was loaded onto a YMC-Triart C18 
basic reversed-phase liquid chromatography column (250 x 4.6 
mm, 5 ym particles) (cat. no. TA12S05-2546WT, YMC, Kyoto, 
Japan). Peptides were separated in a binary buffer system of 2% 
ACN, pH 10 (pH adjusted with ammonium hydroxide, solution 
A), and 98% ACN (without pH adjustment, solution B) in an 
Ultimate-3000 LC system (Thermo Scientific, Massachusetts, 
USA). The gradient of buffer B was set as follows: 0—5% for 1 
min, 5—8% for 5 min, 8—18% for 35 min, 18—32% for 22 min, 
32—95% for 2 min, and 95% for 7 min. Eluates were collected 
every 84 s, and a total number of 50 fractions were obtained. 
The first five fractions were discarded because most of the 
excess labeling reagent was eluted during these fractions. The 
remaining 45 fractions were combined into 15 fractions by 
merging fraction 1, 16, 31; fraction 2, 17, 32; and so on. Then, 
all fractions were dried in a vacuum concentrator and stored at 
—20 °C until further analysis. 


LC-MS/MS Analysis 
For each fraction, an 8 yL sample was injected under trapping 


conditions using an ekspert nanoLC400 system (AB Sciex, 
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Massachusetts, USA) equipped with the ChromXP nanoLC 
trap column (cat. no. 805-00121, AB Sciex, Massachusetts, 
USA) and a ChromXP 3C18-CL-120 column (75 um X 15 cm, 
3 um particles) (cat. no. 805-00120, AB Sciex, Massachusetts, 
USA). Liquid chromatography was interfaced to a tripleTOF 
5600+ instrument (AB Sciex, Massachusetts, USA) via a 
nanoelectrospray ionization source. After being loaded, peptide 
trapping and desalting were carried out at 2 wL/min for 15 min 
with 100% buffer A (2% ACN/0.1% formic acid). The 
analytical separation was established by running a linear 
gradient of buffer B (98% ACN/0.1% formic acid) at a flow 
rate of 300 nL/min from 10 to 27% in 80 min. Then the 
gradient of buffer B was increased to 50% in 24 min and 
decreased to 20% in 6 min. Initial chromatographic conditions 
were restored in 6 s and maintained for 10 min. Data were 
acquired using an ion spray voltage of 2.3 kV, curtain gas of 30, 
gas 1 of 4, gas 2 of 0, and an interface heater temperature of 150 
°C. A rolling collision energy voltage was used for CID 
fragmentation for MS/MS spectra acquisitions. The option 
“adjust collision energy (CE) for iTRAQ” was enabled. Each 
cycle consisted of a TOF MS spectrum acquisition for 250 ms 
(mass range 350—1500 Da), followed by acquisition of up to 30 
MS/MS spectra (100 ms each, mass range 100—1500 Da). The 
total cycle time was fixed to 3.3 s. Dynamic exclusion was set 
for 1/2 of peak width. 


Protein Identification and Quantification 


The LC-MS/MS data for the identification and quantification 
of the global proteome were analyzed using ProteinPilot 
Software 4.5 (AB Sciex, Massachusetts, USA) against the 
Uniprot KB human database (released on June, 2013) with 
88 771 entries, to which 175 commonly observed contaminants 
and all the reverse sequences were added. iTRAQ 8-plex 
labeling and cysteine carbamidomethylation was set as 
modifications. A 95% confidence (unused protein score >1.3) 
was used for protein identification. The PDST_V3.708pE_- 
TEMPLATE (AB Sciex, Massachusetts, USA) was used to 
calculate p values for each ratio. Ratio channels of 114/113, 
116/115, 117/118, and 119/121 were assigned as “target”. 
Ratio channels of 117/119 and 118/121 (technical replicates) 
were assigned as “decoy” (Figure 1). By assigning each ratio 
channel as either “target” or “decoy”, the PDST estimates 
target-decoy global false discovery rate (FDR) for the ratios 
according to their p values. The p value limit was set on the 
basis of target-decoy FDR analysis. Principally, a protein with a 
quantitative ratio differing from 1 and the ratio’s p value less 
than 0.01 is declared as differentially expressed between a pair 
of AS and AI cells. 


Method Validation for MRM Analysis 


Method performance was evaluated for 42 fatty acid methyl 
ester (FAME) standards, eight nucleotides, and acetyl-CoA 
(Supporting Information, Table $1), with respect to linearity 
and reproducibility. For each FAME standard, >0.1 mg/mL 
stock solution was prepared in methanol and stored at —20 °C. 
From these stock solutions, mixtures of 42 FAMEs dissolved in 
hexane were prepared at 0.05, 0.1, 0.5, 1, 5, 10, 50, and 100 ng/ 
mL and analyzed using GC-MS. For each nucleotide and acetyl- 
CoA standard, 1 mM stock solution was prepared in water and 
stored at —20 °C. From these stock solutions, mixtures of the 
nine metabolites dissolved in water were prepared at 0.17, 0.42, 
0.83, 1.67, 8.33, and 83.3 uM and analyzed using LC-MS. 
Linearity was evaluated using the linear regression of peak areas 
with respect to concentration. In addition, FAME standards at 
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Figure 1. Experimental workflow of quantitative proteomics analysis in 
PCa cells. In order to apply the target-decoy global FDR method to 
the protein ratios, two cell lines, Dul45 and 22Rvl, were labeled in 
duplicates. These technical replications function to differentiate the 
real differential detections with false-positive detections, according to a 
p value cutoff. 


0.5 ng/mL were run two times to determine inter-run 
quantitative reproducibility, as measured by the relative 
standard deviation (RSD). 


MRM-Based Relative Quantification of Fatty Acids 


Fatty acids were extracted from PC3, Dul4S, 22Rv1, LNCaP, 
and C42B using a liquid—liquid extraction method.'* Briefly, 
cells were grown to confluence in 6 cm dishes under conditions 
same as those in the proteomic experiment. After the medium 
was removed, cells were washed with 5 mL of cold PBS before 
300 uL of 0.9% KCI solution was added, and cells were then 
detached with a cell scraper. Deuterium-labeled C17:0 fatty acid 
was added as an internal standard. Cells were transferred to an 
Eppendorf Safe-Lock 2.0 mL tube (cat. no. 0030 120.094, 
Eppendorf, Hamburg, Germany). After adding 300 uL of 
methanol and the antioxidant reagent butylated hydroxytoluene 
(0.01% final concentration), cells were incubated overnight at 
—20 °C. Then 1 mL of methyl tert-butyl ether (MTBE) was 
added, and the samples were vortexed thoroughly on a Vortex- 
Genie 2 mixer (cat. no. SI-0236, Scientific Industries, New 
York, USA) for 1 h at room temperature for phase separation. 
After centrifugation at 18 000g for 5 min at 4 °C, the upper 
phase was collected and divided into two parts. One part (200 
HL) was used for analysis of phospholipid-derived fatty acids 
(PFs). The other part was used for analysis of free fatty acids 
(FFs). Both parts were dried under vacuum. PFs were 
derivatizated according to the analytical procedure on the 
Web site of Avanti Polar Lipids company (http://www. 
avantilipids.com/index.php?option=com_content&view= 
article&id=1693&Itemid=412). Samples were dissolved in 1 
mL of hexane. Fifty microliters of sodium methoxide (1 M) was 
added and vortexed thoroughly. Samples were incubated at 
room temperature for 5 min and centrifuged at 1600 rpm at 20 
°C for 5 min. The upper phase was transferred to another Safe- 
Lock tube and vacuum-dried and then dissolved in 50 uL of 
hexane for analysis. For FF derivatization, trimethylsilyl 
(TMS)—diazomethane was used.°° Samples were resolved in 
methyl esterification buffer (500 uL of MTBE/150 pL of 
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MeOH/6.5 uL of 1 M HCl) in a Safe-Lock tube. Fifty 
microliters of TMS—diazomethane (2 M in hexane) was added 
and mixed at a mild speed on a wheel for 30 min at room 
temperature. Six microliters of glacial acetic acid was added to 
quench the reaction. After 1 min, 125 wL of water was added 
and mixed thoroughly. Samples were centrifuged at 1500g for 3 
min for phase separation. The upper phase was collected and 
dried under vacuum. Samples were then dissolved in 50 uL of 
hexane for analysis. Derivatized FFs and PFs were separated 
using Agilent 6890N gas chromatograph (Agilent, California, 
USA) equipped with a J&W HP-88 column (100 m x 0.25 mm 
id. X 0.2 um film thickness) (cat. no. 112-88A7E, Agilent, 
California, USA). The flow of helium carrier gas was set to 1 
mL/min (20 cm/s); the injector temperature was set to 260 
°C, and the temperature of the sample transfer line was set to 
240 °C. Samples were analyzed in splitless mode, and the 
gradient started from 50 °C, held for 1 min, then increased by 
25 °C/min to 175 °C, further by 1 °C/min to 216 and 4 °C/ 
min to 240 °C, then held for 10 min. Mass spectra were 
collected on an Agilent 7000B mass selective detector (Agilent, 
California, USA), using a chemical ionization source in the 
positive mode. Each cell line was analyzed in biological 
triplicate. A washing run was added with every three samples. In 
order to eliminate the influence of fatty acid contaminations, 
three solvent blanks were treated in parallel with samples, and 
the mean values of the peak areas of solvent blanks were 
subtracted from each sample. In addition, a mixture of 24 fatty 
acid standards (Table $3) dissolved in methanol was prepared 
at 5 ng/mL, derivatized with TMS—diazomethane, and used for 
recovery rate evaluation. 


MRM Quantitation of Nucleotides and Acetyl-CoA 


Nucleotides and acetyl-CoA were extracted from five cell lines 
(PC3, Dul45, 22Rvl, LNCaP, and C42B), according to a 
method published earlier.'* Briefly, cells were grown to 
confluence (~10’ cells) in 10 cm dishes under conditions 
that were the same as those in the proteomic experiment, and 
cells were washed twice with ice-cold PBS. Then 100 pL of ice- 
cold water was added into the 10 cm dish, which was then flash 
frozen at —80 °C for about 30 s. Cells were detached with a cell 
scraper on ice and transferred into 1.5 mL Eppendorf tubes. 
Cell suspensions were lysed by two rapid freeze—thaw cycles at 
—80 °C, followed by 30 s of sonication on ice. The volume was 
made up to 400 uL with ice-cold water, and each sample was 
then divided equally into two Eppendorf tubes. Five hundred 
microliters of cold methanol and 250 uL of chloroform were 
added to the homogenates and vortexed. For safety 
consideration, addition of chloroform as well as other toxic 
reagents and sample transfer were always performed in a fume 
hood. Samples were mixed for 10 min at 4 °C and centrifuged 
for 5 min at 2000g at 4 °C to remove debris and precipitate 
protein. The pellet was dissolved in 8 M urea, 0.1 M Tris, pH 
8.5, and protein content was quantified using a BCA kit. The 
supernatant was transferred to a new tube, and 250 pL of 
chloroform and 250 uL of water were added. Samples were 
vortexed for 60 s, kept on ice for 10 min for phase separation, 
and centrifuged for 10 min at 2000g at 4 °C. The resulting 
upper polar and lower nonpolar fractions were separately 
transferred to Eppendorf tubes, vacuum-dried, and combined 
after being dissolved in 30 wL of LC buffer A (5 mM 
ammonium acetate/0.1% FA). A combined volume of a ~60 
uL sample was obtained and centrifuged, and a § mL 
centrifuged sample was analyzed on an Agilent 6460 LC-Triple 
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Quadruple MS (Agilent, California, USA) with an Acclaim 120 
C18 RP column (3 xX 150 mm, 3 um particles) (cat. no. 
063691, DIONEX, California, USA). Metabolites were eluted 
from the column with buffer B (ACN/5 mM ammonium 
acetate/0.1% FA), and buffer A was 5 mM ammonium acetate/ 
0.1% FA. The gradient of buffer B started at 2% and was kept 
for 1 min and then increased from 2 to 30% in 9 min and from 
30 to 80% in 5 min. Then, the gradient of buffer B was 
decreased to 2% in 2 min and kept for 13 min. The total LC 
time was 30 min, with a flow rate of 0.5 mL/min. Gas 
temperature was set at 325 °C, and gas flow was 7 L/min. 
Sheath gas temperature was set at 325 °C, and sheath gas flow 
was 11 L/min. Capillary voltage was 3500 V, and the nozzle 
voltage was 500 V. Delta EMV(+) was 200. Each cell line was 
analyzed in biological triplicate. 


Mitochondrial Staining and Activity Analysis 


Mitochondrial morphology imaging was conducted with 
Mitotracker Green (Life Technologies, Massachusetts, USA) 
as described by Chen et al.'* Cells were seeded in 35 mm Petri 
dishes, and after reaching 50—70% confluence, the cells were 
stained with 50 nM Mitotracker Green diluted in the 
corresponding culture media for 20 min at 37 °C. After 
incubation, cells were washed twice with PBS, stained with 
Hoechst 33342 (Life Technologies, Massachusetts, USA) for 5 
min, and washed three times with PBS. Cell images were 
obtained with an Olympus FV1000 laser scanning confocal 
microscope (OLYMPUS, Tokyo, Japan). 

XF24 Extracellular Flux analyzer (Seahorse Bioscience, 
Massachusetts, USA) was used to detect the extracellular 
acidification rate (ECAR) and oxygen consumption rate 
(OCR) of different cell lines. The medium used for the OCR 
assay was XF assay medium (cat. no. 102352-000, Seahorse 
Biosciences), with 25 mM glucose and 2 mM pyruvate added, 
pH 7.35 + 0.05 (37 °C). The medium used for the ECAR assay 
was DMEM (cat. no. D5030, Sigma), with 143 mM NaCl, 3 mg 
phenol red, and 2 mM L-glutamine added, pH 7.35 + 0.05 (37 
°C). 

Data Analysis 


Enrichment analysis for gene ontology (GO) terms and protein 
interaction network construction was performed using 
STRING. The p values were corrected for multiple testing 
using the method by Benjamini and Hochberg, and the top 10 
significant categories were presented according to their p 
values. Pathway mapping was executed using the Reactome 
database. The Matlab software (version 2010b, MathWorks) 
was used for hierarchical clustering analysis. Proteins were 
clustered using unsupervised hierarchical average linkage 
clustering and a Euclidian distance metric, based on their 
logarithm-transformed ratios. 

For fatty acids, peak areas were integrated with MassHunter 
software (Agilent, California, USA).°° After blank subtraction, 
internal standard correction (this correction was only applied 
for FFs), and protein content correction, the values were log- 
transformed and z-score-normalized and then clustered using 
unsupervised hierarchical average linkage clustering and a 
Euclidian distance metric. Fatty acids with a signal-to-noise 
ratio less than 10 in any sample were excluded for further 
analysis. For nucleotides and acetyl-CoA, the concentrations 
were calculated directly via the standard curves built with 
MassHunter software. After correction for protein content, the 
contents of nucleotides and acetyl-CoA among different 
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Figure 2. Dynamic changes of differentially expressed proteins in different pairs of PCa cells and gene ontology annotations of these proteins. (A) 
Hierarchical clustering based on the 10-base logarithm-transferred ratios of the 528 differentially expressed proteins. By setting the distance cutoff, 
four clusters showing in different colors were identified. (B) Venn diagram depicts the distribution of these differentially expressed proteins in the 
three pairs, and the histogram at the bottom represents the regulation directions of the 191 proteins that change in at least two pairs of PCa cells, for 
which a is PC3-AR*/PC3, b is 22Rv1/Dul45, and c is LNCaP/C42B. (C) Enrichment analysis of gene ontology annotations for the 191 
differentially expressed proteins that changed in at least two pairs of PCa cells, GOBP for biological process, and GOCC for cellular location. 


samples were compared using the student’s t-test with a p value 
of <0.05 as the significance level. 


Safety Considerations 


All procedures involving toxic reagents should be performed in 
a fume hood, and the waste produced must be recovered under 
specific recovery processing. The TMS—diazomethane was 
originally developed as a less toxic substitute for the compound 
diazomethane. However, inhalation of TMS—diazomethane 
may cause lung damage or central nervous system depression 
and should therefore be handled with the appropriate safety 
procedures in place, that is, inside a fume hood, with adequate 
personal safety equipment. Excess TMS—diazomethane is to be 
neutralized using acetic acid. Care must be taken if handled in 
large volumes, as it is highly volatile and produces nitrogen gas. 


M RESULTS AND DISCUSSION 


Differentially Expressed Proteins between AS and Al Cells 


In order to identify the differentially expressed proteins 
between the three pairs of PCa cells, we performed 8-plex 
iTRAQ labeling and used a target-decoy global FDR method to 
set a cutoff p value for the ratios of quantified proteins. A total 
number of 4953 proteins were identified with 1% global FDR at 
the protein level, among which 4847 proteins were quantified 
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(Table S2). When setting the cutoff p value below 0.01, the 
target-decoy global FDR for quantification was 2%. In total, 528 
proteins showed at least one quantitative ratio with a p value 
less than 0.01 (Figure S1). By applying hierarchical clustering 
to the quantitative ratios of the 528 proteins, we found that the 
changes of these differentially expressed proteins were 
heterogeneous among the three pairs of AS and AI cell lines 
(Figure 2A). Among the 528 proteins, 191 were differentially 
expressed in at least two pairs, with 72 proteins increased, 76 
decreased, and 43 proteins changed in both directions in AI 
relative to AS cells (Figure 2B). The other 337 proteins were 
found to be differentially expressed in only one pair (Table S4). 
Only 19 proteins among the differentially expressed proteins 
were changed in the same direction in all of the three pairs, 
among which 10 were increased and 9 were decreased in AI 
relative to AS cells (Figure 2B and Table 1). Two of the 19 
proteins have been shown previously to be critically involved in 
cancer development. They are major vault protein (MVP) and 
poly[ ADP-ribose ]polymerase 1 (PARP-1). MVP was increased, 
and PARP-1 was decreased in AI relative to AS cells. MVP has 
been reported to promote cancer progression by mediating 
epidermal growth factor receptor (EGFR)/PI3K-mediated 
migration and survival pathways in human glioblastoma.” 
Intriguing, abnormal activation of the PI3K pathway is a 
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Table 1. Proteins Showing Consistently Expressional Changes in PCa Cells 


accession protein name 


proteins with increased expressions in AI relative to AS cells 


Q02952-3 isoform 3 of A-kinase anchor protein 12 

P04406 glyceraldehyde-3-phosphate dehydrogenase 
Q9Y2J8 protein-arginine deiminase type-2 

Q14764 major vault protein 

P06737-2 isoform 2 of glycogen phosphorylase, liver form 
P15121 aldose reductase 

Q9BUFS tubulin J-6 chain 

Q14152 eukaryotic translation initiation factor 3 subunit A 
P54577 tyrosine—tRNA ligase, cytoplasmic 

P00338-3 isoform 3 of -lactate dehydrogenase A chain 


proteins with decreased expressions in AI relative to AS cells 


P09874 poly[ADP-ribose] polymerase 1 


P40926 malate dehydrogenase, mitochondrial 

P07954 fumarate hydratase, mitochondrial 

Q86VP6 Cullin-associated NEDD8-dissociated protein 1 

K7EJE8 Lon protease homologue, mitochondrial 

P10515 dihydrolipoyllysine residue acetyltransferase component of 
pyruvate dehydrogenase complex, mitochondrial 

P10809 60 kDa heat shock protein, mitochondrial 

Q06210 glutamine—fructose-6-phosphate aminotransferase 
[isomerizing] 1 

P06744 glucose-6-phosphate isomerase 


common event during PCa progression.’ Here, we also found 
higher expression of EGFR in AI relative to AS cells (Table S4). 
Thus, we speculate that MVP plays a similar role in PCa 
progression as well as in glioblastoma. Previously, PARP-1 was 
found to promote AR transcriptional activity in AR-positive 
prostate cancer cells.'” This could explain the higher expression 
level of PARP-1 in AS relative to AI cells that we observed 
because AS cells exihibit higher AR transcriptional activity 
relative to AI cells. Among the other proteins that were found 
to be changed in the same direction in all three pairs, isoform 3 
of A kinase anchor protein 12 (AKAP12), tyrosine—tRNA 
ligase, cytoplasmic (YARS), isoform 3 of t-lactate dehydrogen- 
ase A chain (LDHA), and glutamine—fructose-6-phosphate 
aminotransferase [isomerizing] 1 (GFPT1) were reported 
earlier as products of androgen-regulated genes,'® in agreement 
with the differences of androgen sensitivity between AS and AI 
cells used in our study. Taken together, these results revealed 
the heterogeneity of PCa cells at the proteomic level. 


Differentially Expressed Proteins Mainly Function in 
Central Metabolism 


To investigate the functions of these differentially expressed 
proteins as well as their roles in CRPC progression, the 191 
proteins whose expressions changed in at least two pairs of AS 
and AI PCa cells were chosen for GO annotation enrichment 
analysis. The result showed that most of these proteins were 
involved in cell metabolism and located in exosomes (Figure 
2C). Exosomes play a pivotal role in cancers,’” *” we thus 
searched the exosome database Exocarta and Vesiclepedia and 
classified these exosome-located proteins based on their 
reported origins (Figure S2). These proteins likely play 
important roles in PCa progression. 

Our results showed that nearly all of the enzymes involved in 
the TCA cycle showed a decrease in expression, while glycolytic 
enzymes including ATP-dependent 6-phosphofructokinase, 
platelet-type (PFKP), glyceraldehyde-3-phosphate dehydrogen- 
ase (GAPDH), a-enolase (ENO1), fructose—bisphosphate 


gene PC3-AR*/PC3 LNCaP/C42B 22Rv1/Dul45 22Rv1/Dul45 
name (114/113) (116/115) (117/118) (119/121) 
AKAP12 0.48 0.08 0.07 0.09 
GAPDH 0.49 0.16 0.54 0.54 
PADI2 0.16 0.36 0.11 0.20 
MVP 0.31 0.22 0.31 0.37 
PYGL 0.42 0.20 0.08 0.10 
AKRIB1 0.19 0.06 0.02 0.03 
TUBB6 0.24 0.14 0.20 0.14 
EIF3A 0.50 0.33 0.67 0.41 
YARS 0.12 0.16 0.41 0.24 
LDHA 0.69 0.08 0.11 0.09 
PARP1 4.33 4.09 7.38 6.43 
MDH2 4.09 10.09 6.37 8.71 
FH 6.61 6.85 5.45 4.79 
CAND1 2.13 2.27 3.02 2.91 
LONP1 5.25 3.66 4.41 5.06 
DLAT 2.21 3.16 2.78 2.13 
HSPD1 1.03 1.15 1.17 1.16 
GFPT1 3.08 3.19 1.80 2.81 
GPI 2.36 1.71 1.56 2 
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aldolase A (ALDOA), and the rate-limiting enzyme LDHA 
showed increased expressions in AI relative to AS cells (Figure 
S3 and Table S4). Specifically, we observed that different 
isozymes of PFK changed in both directions. PFKP was 
increased while PFKM was decreased in AI relative to AS cells 
(Table S4). Previously, it was shown that during neoplastic 
transformation both quantitative increase and isozymic 
alterations of PFK occurred secondary to altered gene 
expression.” In particular, PFK controls the entry of glucose 
metabolites into the downstream steps of glycolysis.” Recently, 
Ros et al. found that the glycolytic enzyme PFKFB4 was 
essential for PCa cell survival,” confirming the importance of 
glycolysis regulation in PCa. 

Key Enzymes Involved in Glutamine Metabolism Change in 
PCa Cells 


Since glucose is mainly utilized during glycolysis, glutamine- 
driven TCA turning powers oxidative phosphorylation and 
accounts for the majority of ATP production in most cancer 
cells.” A recent report showed that PC3 cells use glutamine as 
their main energy source and revealed a key role of glutaminase 
(GLS1) in the growth of PCa cells.” In agreement with these 
studies, we also observed increased levels of GLS1 in C42B 
relative to LNCaP and in Dul45S relative to 22Rvl, and its 
increase in PC3 relative to PC3-AR* was not significant (Table 
S4). However, glutamate dehydrogenase 1 (GLUD1) was 
decreased in C42B relative to LNCaP and in Dul45 relative to 
22Rv1 (Table S4), suggesting a more complex regulation of 
glutamine metabolism in these cells. Recently, GLUD1- 
dependent reductive glutamine metabolism was shown to 
increase in PCa cells when glucose oxidation was decreased by 
metformin.”* The same study observed a more obvious effect of 
metformin in the LNCaP cells, indicating that AS cells were 
dependent more on reductive glutamine metabolism than were 
AI cells. In support of this observation, another key enzyme for 
reductive glutamine metabolism, isocitrate dehydrogenase 
[NADP] cytoplasmic (IDH1), was also shown to be increased 
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Figure 3. Analysis of mitochondrial morphology and function in PCa cells. (A) Mitochondria were stained with MitoTracker Green, and nuclei were 
stained with Hochest 33342. (B) Mitochondria metabolic activity analysis. PC3, Dul45, and C42B cells were plated at 15 000 cells/well, and 22Rv1 
and LNCaP cells were plated at 40 000 cells/well in XF24 V7 cell culture plates 24—48 h prior to the assays. For OCR assay, oligomycin (1 uM), 
FCCP (0.5 uM), and rotenone (1 uM) with antimycin (1 ~M) were added at time points indicated by arrows. For the ECAR assay, glucose (10 
mM), oligomycin (1 yM), and 2-DG (100 mM) were added at time points indicated by arrows. Each data point represents mean + SD (n = 3—4). 
Medium used for different assays was described in the experimental method. For OCR response, data were presented as percentages of baseline 


values (measurement of 3). 


in AS relative to AI cells in our results (Table $4). In addition, 
we found the increase of the mitochondrial aspirate—glutamate 
carrier SLC25A13 in AS relative to AI cells (Table $4). These 
results, together with the information mentioned above, likely 
reflect that glutamine enters multiple metabolic pathways for 
distinct utilities in these PCa cells (Figure SB). 


Functional Changes of Mitochondria Occur between AS 
and Al Cells 


To investigate whether the decrease of enzymes in the TCA 
cycle in AI cells was caused by loss of mitochondrial mass, 
staining of mitochondria and mitochondrial activity analysis 
were conducted in each of the five cell lines (Figure 3). We 
excluded PC3-AR* cells because they showed distinct changes 
of the described TCA enzymes relative to PC3 when compared 
to the other two pairs. We did not observe a significant loss of 
mitochondrial mass in the AI PCa cells from the staining result, 
although the morphology and distribution of mitochondria 


3396 


varied considerably (Figure 3A). Moreover, expressions of the 
mitochondria electron transport chain (ETC) components, 
including electron transfer flavoprotein (ETF) and ATP 
synthase, were increased in AS relative to AI cells except for 
those in the PC3-AR*/PC3 pair in our results (no change of 
ETF and decrease of ATP synthase, Figure $4). However, no 
changes were observed in the other ETC components as well as 
many quantified structural proteins of mitochondria in our 
results, though it may be partially explained by a lack of 
sufficient quantitative peptides. Therefore, functional alter- 
ations rather than mass loss of mitochondria could eventually 
lead to such a complex change of proteins located in this 
organelle among different pairs of PCa cells. To further prove 
this, two mitochondrial proteins showed opposite changes in 
our results. They are the Lon protease homologue mitochon- 
drial LONP1 that increased and calcium-binding mitochondrial 
carrier protein SCaMC-1 (SLC25A24) that decreased in AS 
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Figure 4. Differences in the contents of nucleotides and acetyl-CoA in PCa cells. The cellular contents of ADP, ATP, AMP, and cAMP in PCa cells 
are shown as mean + SD (n = 3) in (a), and their ratios to ATP are shown as mean + SD (n = 3) in (b). The cellular content of acetyl-CoA is shown 
as mean + SD (n = 3) in (c). The cellular content of NAD and NADH is shown as mean + SD (n = 3) in (d), and the cellular content of NADP and 
NADPH is shown as mean + SD (n = 3) in (e). The ratios of NADH/NAD, as well as NADPH/NADP, are shown as mean + SD (n = 3) in (f). The 


letters above the columns denote significant differences (p < 0.05) for comparison, among which 


LNCaP and “b” represents comparison between Dul45 and 22Rv1. 


relative to AI cells (Table S4). LONP1 plays a key role in 
metabolic reprogramming by remodeling OXPHOS com- 
plexes.” SLC25A24 regulates the contents of adenosine 
nucleotides in the mitochondrial matrix and increases the 
antiapoptotic activity of mitochondria.” Both of them are 
involved in timely response of mitochondria to cellular 
metabolic dynamics. An even more complex change was 
observed for dynamin-like 120 kDa protein, mitochondrial 
OPAI, which was increased in PC3-AR* relative to PC3 and in 
LNCaP relative to C42B, but was decreased in 22Rv1 relative 
to Dul45 (Table $4). OPA1 is the key effector for fusion of the 
mitochondrial inner membrane. Importantly, the rate of 
mitochondrial fusion responds to changes in metabolism.* 

Together, these results reveal that metabolism-related 
functional alteration rather than mass loss of mitochondria 
occurs between AS and AI cells. 


High Rate of Glycolysis Serves the High Demand for a 
Biosynthetic Precursor in Al Cells 


Cancer cells are known to rewire their metabolism and energy 
production networks to support rapid proliferation.” In 
contrast to resting cells, proliferating cancer cells are in much 
greater need of building blocks for reductive biosynthetic 
reactions.*>°> 

To further investigate the metabolic signatures of PCa cells, 
we compared the contents of several metabolites in PCa cells 
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“n 


a” represents comparison between C42B and 


by the MRM approach. In Table S1, we summarize the results 
of MRM method vaidation for all targeted metabolites in this 
study. All metabolites gave a linear response (R° > 0.97). As 
shown in Figure 4 and Table S5, the contents of ADP, AMP, 
cAMP, NAD, NADH, and NADP, as well as the ratios of ADP/ 
ATP, AMP/ATP, cAMP/ATP, NADH/NAD, and NADPH/ 
NADP, showed significant differences in one or two pair of PCa 
cells. Since the decrease of TCA enzymes in AI relative to AS 
cells was observed in our proteomic results, the increased 
content of NADH and the increased ratio of NADH/NAD in 
AI relative to AS cells might be caused by other mechanisms. 
One possible mechanism is the regeneration of NAD via the 
lactate generation step in glycolysis, as a high rate of glycolysis 
marks AI cells (Figure 3B). The malate—aspartate shuttle 
(Figure 5B) is another pathway that regenerates NAD but is 
kinetically much slower than pyruvate-to-lactate conversion.” 
Therefore, elevated glycolysis will lead to a higher ratio of 
NADH/NAD in AI relative to AS cells. 

Since the activity of PFK was reported to be sensitive to the 
AMP/ATP ratio,* and we also observed a higher ratio of 
AMP/ATP in AS relative to AI cells, PFK could control the rate 
of glycolysis via isoenzymatic regulation or activity change 
according to AMP/ATP ratio in PCa cells. 

More importantly, the elevation of the AMP/ATP ratio as 
well as the increased content of AMP in AS cells could lead to 
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Figure 5. Relative comparison of fatty acids and the changes of related metabolic enzymes in PCa cells. (A) Heatmaps of the free fatty acids (FF) 
and phospholipid-derived fatty acids (PF) detected in different PCa cell samples; the numbers following each cell line’s name represent biological 
replicates. The contents of short-chain PFs (less than 20 carbons), as well as polyunsaturated C20:3n6 and C20:4n6 PFs, were higher in AS relative 
to AI cells (the green cluster), similar to the changes of FFs (the blue cluster), explaining them clustering together. The contents of long-chain PFs 
were decreased in AS relative to AI cells (the orange and red clusters). The most consistent changes were observed for several monounsaturated PFs, 
namely, C18:1n9c, C20:1n9, C22:1n9, and C24:1 (the red cluster), which were increased in AI relative to AS cells. Polyunsaturated C22:5n3 and 
C22:6n3 PFs were specifically upregulated in C42B cells, similar to the change of C22:0 PF (the orange cluster). (B) Schematic cellular metabolic 
pathways and the altered enzymes shown in colors in PCa. The green dashed box highlights fatty acid synthesis and beta oxidation that serve as the 
major energy supply in AS cells. MUFA, monounsaturated fatty acid; PUFA, polyunsaturated fatty acid; PLFA, phospholipid fatty acid; MAG, 


monoacylglycerol; ChE, cholesteryl ester. 


AMPK activation. AMPK is a core regulator of cellular 
metabolism. A previous study using human PCa specimens 
showed that AMPK is highly activated in 40% of the cases, and 
it serves as a downstream effector of AR-induced PCa cell 
proliferation.*© In androgen-treated PCa cells, the androgen- 
mediated AMPK-PGC-1a signaling cascade is responsible not 
only for increasing the rate of glycolysis but also for increasing 
glucose and fatty acid oxidation.” Furthermore, Ca”*/calm- 
odulin-dependent protein kinase kinase # (CaMKK/), an 
AMPK kinase, was found to be highly expressed in normal 
prostate tissue and elevated in prostate cancer.’ One recent 
study reported that a novel AMPK activator inhibited the 
growth of PCa cells.’ It is worth mentioning that the study 
includes different AR phenotypes, and the basal level of AMPK 
activation and the functions of AMPK activator are different in 
these PCa cells. However, the role of AMPK in prostate cancer 
is controversial.*°*! The timing, level, and cellular context of 
AMPK activation may have profound effects on the functional 
consequences of AMPK signaling.” A more complicated 
negative feedback loop exists between AR and AMPK in PCa 
cells} 

In our proteomic results, an increased expression of AK1 in 
AS relative to AI cells was observed (Table S4). The highly 
reversible adenylate kinase reaction can be partially displaced 
toward synthesis of AMP and cause a large rise in AMP,** 
representing an alternative route for AMPK activation in AS 
cells (Figure 6). Further studies are required to confirm these 
results. 
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Simultaneous Increase of Fatty Acid Synthesis and Fatty 
Acid Oxidation in AS Cells 


A metabolic shift, in which PCa cells use alternative enzymes 
and pathways to facilitate the production of fatty acids, is 
proposed to contribute to the well-established lipogenic 
phenotype of PCa.** In agreement with this, we observed an 
increase of fatty acid synthase (FASN) in AS relative to AI cells 
(Table S4). Importantly, FASN has been reported to be a 
candidate oncogene in PCa.“ However, changes of other 
enzymes in fatty acid synthesis, such as ATP-citrate synthase 
(ACLY) and acetyl-CoA carboxylase 1 (ACC), were complex 
among the different pairs (Table $4). ACLY was decreased in 
LNCaP relative to C42B but increased in 22Rv1 relative to 
Dul45. ACC was only increased in 22Rv1 relative to Dul45 
and showed no difference in the other two pairs. In agreement 
with the changes observed for ACLY and ACC, the content of 
acetyl-CoA was higher in C42B cells compared to LNCaP as 
well as to those in the other PCa cells, but showed no 
difference between Dul45 and 22Rv1 (Figure 4c), suggesting 
that acetyl-CoA accumulated in C42B but was consumed in 
22Rv1 through fatty acid synthesis. Another function of acetyl- 
CoA is to serve as substrate for cholesterol synthesis. We 
observed that squalene synthase (FDFT1), the key enzyme for 
cholesterol synthesis, was present at an elevated level in AS 
relative to AI cells (Table S4), a fact that is in agreement with 
the increased de novo synthesis of cholesterol in PCa cells.*” 
Recently, AI cells were shown to bear aberrant accumulation of 
esterified cholesterol in lipid droplets, which arises from an 
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Figure 6. Proposed model of the AR signaling network and its role in cellular metabolism regulation. Different status of AR activity regulates the 
transcription of a distinct set of androgen-responsive genes (ARGs) and results in the two types of PCa cells, that is, AS and AI (also called CRPC 
cells). A major role of AR is its inhibitory effect on in situ androgen synthesis, which decreases intracellular androgen level directly and forms a 
negative feedback loop in AR activity regulation. This negative feedback loop suppresses AR signaling at high androgen levels but allows increased 
AR and androgen synthesis in CRPC. ADT induces the change of intracellular androgen level, switching AR to partial activation, and finally 
contributes to androgen independence. AR activation will increase the expression of CaMKKf and activate AMPK, and the latter will conversely 
inhibit AR activation. As a central regulator of cellular metabolism, AMPK mediates metabolic shift in PCa cells. Different ARGs regulated by distinct 
AR activity status also contribute to this shift. AK1 can convert ADP to AMP and increase the AMP/ATP ratio, leading to AMPK activation. 
Different metabolic pathways highlighted in the colored octagonal boxes are dominant in distinct PCa cells, and the related enzymes that are 
differentially expressed in PCa cells are highlighted in the colored boxes above each pathway. 


uptake of exogenous lipoproteins and is independent of 
androgen signaling.“ Thus, AI cells appear to obtain 
cholesterol through exogenous uptake rather than de novo 
synthesis, which is in contrast to AS cells. 

After synthesis, fatty acids can be oxidized for energy supply 
or converted into the carbon skeleton of triglyceride for lipid 
storage as well as for phospholipid synthesis.” We observed 
that carnitine O-palmitoyltransferase (CPT) changed in the 
isoform context among different pairs of PCa cells (Table $4). 
CPT1A and peroxisomal carnitine O-octanoyltransferase 
(CROT), which are located on the outer membrane of 
mitochondria and peroxisomes, respectively, were found to be 
increased in LNCaP relative to C42B, while showing no 
difference in the other pairs. CPT2, which locates on the inner 
membrane of mitochondria, was increased in LNCaP relative to 
C42B and in 22Rv1 relative to Dul4S. One reason for this 
observation could be that oxidation for different fatty acids in 
mammalian cells is compartmentalized according to the lengths 
of their acyl chains. In comparison with AI cells, AS cells 
overexpressed several mitochondrial enzymes that are respon- 
sible for fatty acid J-oxidation, including long-chain fatty acid— 
CoA ligase 3 (ACSL3), enoyl-CoA hydratase (ECHS1), 
Delta(3,5)-Delta(2,4)-dienoyl-CoA isomerase, mitochondrial 
(ECH1), trifunctional enzyme subunit alpha (HADHA), and 
acetyl-CoA acetyltransferase (ACAT1) (Table S4). Consistent 
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with the increased levels of these enzymes, most of the 17 FFs 
reliably detected by MRM were found to be increased in AS 
relative to AI cells (Figure SA, blue cluster). Together, these 
results indicate that fatty acid -oxidation in mitochondria is 
increased in AS cells. This may result from the androgen- 
mediated AMPK-PGC-la@ signaling cascade, which directly 
causes the increase of fatty acid f-oxidation.*” In order to 
investigate the fates of fatty acids in PCa cells, we also 
compared the changes of PFs in different PCa cells relatively 
using MRM. Consistent with the complex changes of lipogenic 
enzymes observed here (data not shown), the contents of 
different PFs showed complex changes between AS and AI cells 
(Figure 5A). Our results were consistent with a recent MS 
study on phospholipid composition in PCa cells.°° That study 
reported that AS cells contain more polyunsaturated 
phospholipids but less saturated and monounsaturated 
phospholipids than AI cells. Moreover, we observed that the 
contents of short-chain PFs as well as FFs were higher in AS 
relative to AI cells, indicating that synthesized FFs were 
partially converted to short-chain PFs in AS cells. Long-chain 
PFs, especially the monounsaturated counterparts, displayed 
changes in opposing directions in our results, suggesting 
preferential uptake from the environment rather than de novo 
synthesis in AI PCa cells. 
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Thus, we suggest that a simultaneous increase in synthesis 
and /-oxidation of fatty acids occurs in AS cells, which may 
serve as a major energy supply (Figure 5B). In agreement with 
our hypothesis, a previous study showed that, in PCa cells, 
androgen could induce the activation of both lipogenesis and 
fatty acid degradation in a coordinated fashion.°’ Further 
studies will be required to clarify the regulation of lipid 
metabolism in CRPC cells. 


Al PCa Cells Favor Cell—ECM Focal Adhesion over Cell—Cell 
Adhesion 


Another category of proteins that was observed to show 
changes in expression between AS and AI is that of factors 
associated with cell adhesion. These factors included tubulin /- 
6, talin-1, Fermitin family homologue 2, IQGAP1, fascin, actin- 
related protein 3, CAP1, TWF2, coronin, and profilin-1. All of 
these proteins were found to be present at elevated levels in AI 
relative to AS cells. In addition, @-actinin-1 and isoform 1 of 
vinculin changed in the same direction between AS and AI cells 
(Table $4). AI cells also expressed more integrin 1 and less E- 
cadherin relative to AS PCa cells (Table S4). It is well- 
established that integrin #1 mainly mediates cell adhesion to 
extracellular matrix (ECM) components, whereas E-cadherin 
mainly mediates epithelial cell-cell adhesion.°* According to 
the model of adhesion dynamics proposed by Parsons et al.,>* 
we conclude that AI PCa cells favor focal adhesion to ECM 
over cell—cell adhesion based on the changes of these adhesion- 
related molecules observed here. Adhesion-mediated signaling 
through integrin #1 and integrin-linked kinase has a critical role 
in caveolae formation.** In addition to caveolae formation, 
caveolin-1 plays many other roles in CRPC progression.'° 
Consistently, we observed that caveolin-1 was increased in AI 
relative to AS cells, though the p values were above 0.01 (data 
not shown), which may be a result of insufficient numbers of 
peptides for quantification. Another key cell-ECM adhesion 
molecule, CD44, was also increased in AI relative to AS cells 
(Table S4). Previously, CD44 was shown to connect to the 
actin cytoskeleton, bind to growth factors, and enhance 
signaling through cis interactions with EGER.” Whether 
CD44 links to EGFR for PCa cell adhesion needs further 
investigation. 

Overall, AI cells favor integrin-/1-mediated cell adhesion to 
ECM over E-cadherin-mediated cell—cell adhesion. Like the 
metabolic enzymes, most of these adhesion-related molecules 
also locate in the exosome (Figure $2). Therefore, elucidating 
the internal relations between cell adhesion and metabolism 
may finally uncover the role of exosomes in PCa. 


m CONCLUSION 


An increasing number of studies reveal the heterogeneity of 
PCa as a problem for the identification of common features in 
its biological properties. In this study, we took this 
heterogeneity of PCa into consideration and compared the 
proteomes of three pairs of PCa cells. We identified 19 
differentially expressed proteins that changed in the same 
direction in all three pairs, among which four proteins have 
been reported to be the products of androgen-responsive genes. 
Our results suggest that the 19 proteins are pivotal in CRPC 
progression. Furthermore, isozymic alteration of PFK, together 
with the resulting aberrant glycolysis, may accompany CRPC 
progression, based on the established role of glycolysis in rapid 
proliferating cancer cells.” Apart from altered glycolysis, 
glutamine and fatty acid metabolisms are also different between 
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AS and AI cells according to our proteomic results as well as 
MRM analysis of fatty acids, nucleotides, and acetyl-CoA. 
These changes suggest a metabolic shift between AS and AI 
cells and reveal a central place of mitochondria in this shift. 
Increased glycogen mobilization is likely to represent an 
additional sign for the metabolic shift, based on our 
observations that two isoforms of glycogen phosphorylase, 
namely, PYGL and PYGB, were elevated in AI relative to AS 
cells (Table 1 and Table S4). The simultaneous increases of 
AMP and AK] likely contribute to AMPK activation in AS cells, 
suggesting a potential mechanism underlying the metabolic 
shift (Figure 6). Furthermore, integrin-f1-mediating cell 
adhesion to ECM was preferentially utilized in AI relative to 
AS cells. Because many differentially expressed metabolic 
enzymes as well as adhesion molecules observed in our study 
are located in exosomes, the synergistic action between cell 
adhesion and metabolism in exosomes confirmatively contrib- 
utes to CRPC progression and deserves further investigations. 
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